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ABSTRACT 

Cosmic dust extinction alters the flux of type la supernovae. Inhomogeneities in the dust distribution 
induce correlated fluctuations of the SN fluxes. We find that such correlation can be up to 60% 
of the signal caused by gravitational lensing magnification, with an opposite sign. Therefore if not 
corrected, cosmic dust extinction is the dominant source of systematic uncertainty for future SNe la 
lensing measurement limiting the overall S/N to be < 10. On the other hand, SN flux correlation 
measurements can be used in combination with other lensing data to infer the level of dust extinction. 
This will provide a viable method to eliminate gray dust contamination from the SN la Hubble 
diagram. 

Subject headings: Cosmology: large scale structure-gravitational lensing- supernovae: general-dust: 
extinction 
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1. INTRODUCTION 

Gravitational lensing causes several observable effects 
such as distortion of galaxy shape (cosmic shear), vari- 
ation of galaxy number density (cosmic magnification) 
and mode-coupling in cosmic backgrounds. Over the 
upcoming years measurements of these effects will pro- 
vide an accurate mapping of t he matter distribution in 
the u n iverse (for revie ws see iBartelmann fc Schneider! 
(|Mnh : lRefregIer1 (pOOlh L 

Recently, several other lensing reconstruction methods 
have been proposed. One possibility is to measure the 
spatial correlation of lensing induced supernova (SN) flux 
fluctuations. In fact due to lensing magnification 1 , the 
SN flux is altered such that F — ► F^i ~ F(l + 2k), where 
F is the intrinsic SN flux, fi is the lensing magnifica- 
tion and k is the lensing convergence. Intrinsic fluctu- 
ations of the SN flux are random (analogous to intrin- 
sic galaxy ellipticities in cosmic shear measurement). In 
contrast those induced b y lensing magnifi c ation (see e.g . 
i Kantowski et al. I (| 1995( 1 : iFrieman I (|1996h : iHoizI ([1998( 1; 
I Dalai et alj (|2003D ) are correlated with the overall mat- 
ter distribution (analogous to the shear signal). There- 
fore the lensing signature can be inferred ei ther from spa- 
tial c orrelation measurements of SN fluxes ([Cooray et al.l 
2006) or from the root-mean-square of flux fluctuations 
of hig h redshift SNe for which the l ensing signal is dom- 
inant (|Dodelson fc Vallmottoll2006f) . 

Gravitational lensing also induces scatter in the galaxy 
fundamental plane through magnification of the effective 
radius, R e — » it^/i 1 / 2 ~ i? e (l + «)■ Since intrinsic scat- 
ters in the fundamental plane are random, spatial cor- 
relatio n measurements can be u sed to infer the lensing 
signal (Bertin & Lombardi 2006). A similar analysis can 
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1 Throughout this paper, the term lensing magnification refers 
to both the cases of magnification (/i > 1) and de-magnification 
(fi < 1). To be more specific, the spatial correlation functions 
and the corresponding power spectra (C K and C K s T ) investigated 
hereafter are averaged over the full distribution of fi. 



be applied to the Tully-Fisher relation as well. 

Astrophysical effects may limit the accuracy of these 
methods. For instance extinction by cosmic gray dust 
can be an important source of systematic uncertainty. 
This is because dust absorption changes the apparent 
SN flux and may induce correlation of the flux fluctua- 
tions. It also induces scatters in the fundamental plane 
by dimming the galaxy surface brightness and affects the 
Tully-Fisher relation through dimming the galaxy flux. 
These effects potentially cause non-negligible systematics 
in the corresponding lensing measurements. 

Although the existence of gray dust in the inter- 
galactic medium (IGM) remains untested, this scenario 
could account for the m etal enrichment of the IGM 
([Bianchi fc Ferraral (|2005T) and reference therein). Test- 
ing the gray dust hypothesis is also relevant for cosmolog- 
ical parameter inference from SN la luminosity distance 
measurements. Recently ICorasanitil (j2006[ ) has pointed 
out that gray dust models which pass current astrophys- 
ical constraints can induce a ~ 20% bias in the estimate 
of the dark energy equation of state w using the Hubble 
diagram of future SN la experiments. 

In this paper, we study the impact of cosmic gray dust 
on SN lensing measurements, under the optimistic as- 
sumption that contaminations of reddening dust can be 
perfectly corrected. The effects on lensing reconstruc- 
tion based on the fundamental plane and the Tully- 
Fisher relation can be estimated similarly. For super- 
nova, the key point is that extinction caused by dust 
inhomogeneities along the line of sight causes flux fluc- 
tuations which are anti-correlated with the lensing mag- 
nification signal and thus wash-out its imprint. In par- 
ticular we find that dust induced correlation can bias SN 
lensing measurements by 10 — 60%. Therefore this ef- 
fect is likely to be the dominant source of systematics 
for future SN surveys characterized by large sky cover- 
age and sufficiently high surface number density. If not 
corrected, the dust induced correlation would limit the 
signal-to-noise to S/N < 10. This is low compared to 
the S /N achieved by current cosmic shear measurements 



("e.g. Uarvis et al I (f2005h ; IVan Waerbeke et all (|2005l ); 
iHoekstra et al.l ( 2005T )) an d that of proposed meth- 
ods such as CMB lensing (ISeliak fc Zaldarriaga I Il999t 
Zaldarriaga fe Seliak I [1999; Hu 1 120011: jHu fe Okamotol 
2002ft. 21cm background lensing (ICooravll2004HPerJl2004 
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and cosmic magnification of 21cm emitting galaxies 
(|Zhang fc Pen 1120051 l2O0l . 

Nevertheless we suggest that measurements of the SN 
flux correlation still carry valuable information. In fact 
in combination with other lensing data they will pro- 
vide a viable method to detect and eliminate cosmic 
gray dust contamination from future SN la luminosity 
distance measurements. 

2. DUST INDUCED FLUX FLUCTUATIONS 

The observed flux of a SN la at redshift z in the direc- 
tion n of the sky is 

7obs 



F obs (h,z) = Fjie 



(1) 



where F is the intrinsic flux and r the optical depth 
caused by dust extinction along the line of sight. The 
lensing magnification can be written as /i = 1/[(1 — k) 2 — 
7 2 ] ~ 1 + 2k with k and 7 being the lensing convergence 
and shear respectively. In the presence of dust density 
inhomogeneities, the optical depth can be decomposed in 
a homogeneous and isotropic part f and a fluctuation 5t 
(t = f + St). To first order, Eq. ([1]) reads as 

F obs (n, z) ~ Fe~ f {z] [1 + 2n(n, z) - 5t{ti, z)} . (2) 

It is worth noticing that the lensing and dust extinction 
terms have opposite sign. Since /2 = (ensemble average) 
and R = 0, the average flux of a SN la sample in a given 
redshift bin is F ohs (z) ~ Fe~ f ^\ 

The angular correlation of the flux fluctuations can be 
inferred from the e stimator 5p{n 1 



— z^obs 



pons I jpobs ^ 

(jCoorav et al.ll2006f ). From Eq. ^ we then have S F = 
2k, — St, hence Sf provides an estimate of the gravita- 
tional lensing only if fluctuations in the optical depth are 
negligible. 

The lensing convergence k is related to the 3D matter 
overdensity 5 m by 



S m W(x,Xs)dx , 



(3) 



where W(\tXs) is the lensing geometry function. For a 
flat universe W(x> Xs) = (1 + - z )x(l — x/Xs), with x an d 
Xs the comoving diameter distance to the lens and source 
respectively. 

Following the derivation of lCorasanitil (|2006f ) the aver- 
age optical depth to redshift z is 



f(z) 



1 



2.5 log e 
with c the speed of light and 



dA 
dz 1 



cdz , 



(4) 



1 



dA 



Pd(z) 



2.51ogedz Ag(l + z)H(z) 



Qi(a,z)N(a)—, 
a 



(5) 

where pd is the average dust density, g is the grain ma- 
terial density, a is the grain size, is the extinction 
efficiency factor at the rest-frame wavelength A which 
depends on the grain size and complex refractive index 
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Fig. 1. — The lensing normalized matter surface density and 
the B-band dust extinction Ag for different dust models (see text). 
Since Ab and are comparable, dust extinction effects cannot 
be neglected in lensing measurements of SN flux correlation. 



m, and N(a) is the size distribution of dust particles. 
The extinction efficiency factor is computed by solving 
numerically the Mie equations for spherical grains (Bar- 
ber & Hill 1990). Since dust particles are made of met- 
als we estimate the evolution of the average cosmic dust 
density pd from the redshift dependence of the average 
cosmic metallicity as inferred by integrating the star for- 
mation history (SFH) the Universe. Such a modeling 
is an extention of that presented in Aguirre (1999) and 
Aguirre & Haiman (2000), since in addition to estimat- 
ing the amount of cosmic dust density in terms of the 
measured SFH, it accounts for the physical and optical 
properties of the dust grains. 

This approach differs from that used in some of the 
SN la literature (see for instance Riess et al. (2004)). 
In these studies the cosmic dust dimming is estimated 
by modeling the evolution of dust density as a redshift 
power law with different slopes corresponding to different 
cosmic dust models. More importantly these studies as- 
sume the empirical interstellar extinction law, tipically in 
the form inferred by Cardelli et al. (1989). However cos- 
mic dust particles undergo very different selection mech- 
anisms compared to interstellar grains and therefore are 
unlikely to cause a similar extinction. 

In this perspective our modeling is rather robust, since 
the cosmic dust absorption is computed from first princi- 
ples and in terms of astrophysical parameters which can 
be measured through several observations, such as X-ray 
quasar halo scattering (see Paerels et al. 2002) or high 
resolution measurements of the Far Infrared Background 
(FIRB) (Aguirre & Haiman 2000). For more details on 
these cosm ic dust models an d their cosmological impact 
we refer to ICorasanitil (|2006l ). 
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The fluctuation in the optical depth is then given by 



(6) 



where Sd is the fractional dust density perturbation. The 
resulting auto-correlation power spectrum of Sf is: 



(7) 



where C K , C$ T , C K $ T are the angular power spectra of 
k, St, and the k-5t cross correla tion. Using the Lim- 
ber's approximation these read as ([Limber Ifl9 54: Kaiser 
1998): 



l 2 C K = n 
2tt I 

l 2 C Sr 7T 



2tt 



2c 2 



1 

2.5 log e 



J A 2 (~,z\ W 2 ( X ,Xs)xdx 



(8) 



dA 



Xdx 



and 



(9) 



l 2 C K s T 7T 3Q m H 2 

I 5c 2 log e 



2tt 



A 2 
^ss d 



-,z] W(x,Xs)^r-xdX, 
X J ax 

(10) 

where A 2 = k 3 Ps(k)/2w 2 is the dimcnsionlcss matter 
density variance and P$ is the matter density power spec- 
trum. The nonlinear A 2 is calculated using the Peacock- 
Dodds fitting formula (jPeacock k Doddd 119961 ). A 2 SSd 
and A 2 ^ are defined analogously. The spatial distribu- 
tion of IGM dust is not known, the simplest assumption 
is that dust traces the total mass distribution. In such 
case A 2 d = b^A 2 and A 2 Sd = b^A 2 , where bd is the dust 
bias. 

Defining T, L = f^mff fW{x_,Xs)dx, one has St/k ~ 
b d A/X L , hence C 5r /C R ~ b 2 {A/^ L ) 2 and C k5t /C k - 
bd(A/YjL). This indicates that cosmic dust contamina- 
tion is negligible only if A(z) <C Sl(z). 

We adopt a flat ACDM cosmology with fl m = 0.3, 
fi A = 0.7, h = 0.7, Q b = 0.04, cr 8 = 0.9 and the primor- 
dial powe r index n = 1, We a ssume the BBKS transfer 
function (|Bardeen et al.lH986l ). For the dust extinction 
we limit our anal ysis to a te s t-bed of four cosmic dust 
models studied in ICorasanitil (|2006f ). These are charac- 
terized by model parameter values motivated by astro- 
physical considerations. In particular the particle size 
distribution is in the range 0.05 — 0.2 fim, consistently 
with the fact that smaller grains are distroyed by sput- 
tering, while larger ones remain trapped in the gravita- 
tional potential of the host galaxy (Ferrara et al. (1991); 
Shustov & Vibe (1995); Davis et al. (1998)). The grain 
composition consisting of Silicate or Graphite particles 
and we consider both low and high star formation his- 
tory scenarios. 

Models A and B refer to Graphite particles with 
low and high SFH respectively, while models C and D 
correspond to Silicate grains. The total dust density 
for these models is within the limits imposed by the 
DIRBE/FIRAS data (Aguirre & Haiman 2000) and con- 
cides with the upper limit obtained from the analysis of 
X-ray quasar halo scattering (Paerels et al. 2002). These 
gray dust models cause little reddening of the incoming 
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Fig. 2. — Lensing and dust contamination power spectra. The 
upper line is the lensing convergence 1 2 C k /2tt. Other lines are 
CkSt ~C$t/^ f° r dust model A, B, C and D, respectively with = 
1. We have assumed all SNe to be at z s = 1. Clearly, the existence 
of cosmic dust would degrade or even prohibit measurement of the 
lensing signal. 



light and induce a color excess in the optical and near-IR 
bands smaller than 0.01 mag. 

A further assumption concerns the gray dust spatial 
distribution, for which we have little knowledge of. This 
model uncertainty may affect the results presented in this 
paper significantly. One can imagine an extreme case 
where gray dust distributes homogeneously. Then there 
will be no fluctuations in r and thus no induced corre- 
lation in SN flux fluctuations. However, current under- 
standing of gray dust formation implies that gray dust 
is associated with the overall matter distribution. So 
a more appropriate treatment of gray dust distribution 
is the bias model 8^ — bd8 m , as adopted in this paper. 
Although it is natural to expect bd to be redshift and 
scale dependent, since we have little knowledge of it for 
simplicity we assume bd = 1. 

From Fig. Q] we can see that S l is comparable to the 
B-band dust extinction A, hence dust contamination can- 
not be neglected. Therefore C$ T and C K $ T in Eq. ([7]) are 
source of systematic errors which need to be corrected if 
we aim to measure the convergence power spectrum. 

In Fig. [2] we plot the lensing convergence power spec- 
trum l 2 C K /2ir, and the dust contamination power spec- 
trum C k St — C(5t/4 for our test-bed of dust models for 
sources at z s = 1. We find that Cg T is smaller than C K $ T , 
mainly due to the 1/4 prefactor. Since C K g T has an op- 
posite sign to the lensing signal in Eq. (O, its overall 
effect is to suppress the spatial correlation of SN la flux 
fluctuations and consequently diminish the variance and 
covariance of flux fluctuations. Since statistical errors on 
cosmological parameters constraints from SNe la Hubble 
diagram are proportional to the square root of the vari- 
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Table 1 

The relative error caused by extinction with 

RESPECT TO LENSING, Tj = \C K Sr — Cs T /4\/C K AT 
DIFFERENT SOURCE REDSHIFTS FOR OUR TEST-BED OF DUST 
MODELS. 77 IS ROUGHLY INDEPENDENT OF MULTIPOLE I, 
SINCE SHAPES OF C K , C K g T AND Cj T ARE VERY SIMILAR. 



source redshift 

z s 


Graphite 
high SFH 
Model A 


Graphite 
low SFH 
Model B 


Silicate 
high SFH 
Model C 


Silicate 
low SFH 
Model D 


0.5 


0.65 


0.33 


0.24 


0.11 


1.0 


0.45 


0.20 


0.21 


0.08 


1.7 


0.36 


0.13 


0.19 


0.06 



ance and covariance (see e.g. ICoorav et al.l (|2006af) for 
discussions), the existence of cosmic dust extinction fluc- 
tuations decreases the statistical uncertainties, though 
the mean dust extinction will induce a systematic bias 
unless corrected. 

Dust contamination can be quantified by the ratio 
77 = \C K $ T — Cs T /4:\/C K . Since both k and St trace the 
same large scale structure (enforced by the simplification 
bd =const.), the multipole dependence of C K , Cs T and 
C K $ T are similar such that 77 is roughly constant. In ta- 
ble[T]we list its values for sources at redshift z s — 0.5, 1.0 
and 1.7 respectively. As it can be seen, model A causes 
the largest contamination inducing a systematic error as 
large as 60% of lensing signal. Even for model D the 
contamination is still ~ 10%, which is comparable to the 
statistical error expected from future SN la lensing mea- 
surements. Consequently dust induced systematics will 
be the dominant source of uncertainty for this type of 
measurements . 

Furthermore we find that the relative error can be ap- 
proximated by 77 = (3bdA/Y,L, where (3 ~ 0.7 with a 
dispersion < 0.1 over the redshifts investigated for our 
test-bed of dust models. This relation suggests that if we 
can measure 77 in combination with an independent lens- 
ing measurement, it would be possible to infer A given 
knowledge of bd- In the next section we will discuss how 
these type of measurements can be used to remove cos- 
mic dust contamination in the SN la Hubble diagram. 

3. REMOVING COSMIC DUST CONTAMINATION 

Flux fluctuations induced by lensing and extinction 
are small compared to intrinsic SN flux fluctuations and 
therefore can only be extracted statistically, except for 
the strongly lensed or heavily extincted SNe. Accurate 
lensing measurements can be obtained from a variety of 
astrophysical observations of cosmic shear and cosmic 
magnification. In combination with correlation measure- 
ments of SN fluxes these can be used to quantify the level 
of cosmic dust extinction and provide a viable method to 
remove dust systematics from the SN la Hubble diagram. 
The idea is to infer rj from the comparison of Cg F and 
C K . As discussed_ before 77 ~ Q.7bdA/T,L, this would al- 
low to measure A up to model uncertainties in bd and 
measurement errors in Cs F ■ The estimated value of A 
can then be used to correct the standard-candle relation 
of SN la. 

The efficiency of this method depends on the sky cov- 
erage and the SN number density of the survey. For 



instance in order to measure A to 10% accuracy, the 
overall S/N of C Sp must be > 10(1/7? - 1). This im- 
plies that for model A, C$ p should be measured with 
S/N of ~ 10, while for model B, C and D, it would re- 
quire a S/N > 40-100. In the case of a survey with 10 4 
SNe and covering 20 deg 2 of the sky the signal-to-noise 
is S/N ~ 10 (ICoorav et al.ll2006h . Since S/N oc f)£ 
to reach S/N — 40-100 requires a factor of 20-100 times 
higher in sky coverage and total number of observed SNe. 
This c an be achieved by the proposed ALPACA experi- 
ment (|Corasaniti et al.ll2006fh 

Galaxy-quasar correlation measurements provide an- 
other method to estimate the level of cosmic dust ex- 
tinction. For a given line of sight, dust extinction re- 
duces the observed number of galaxies above flux F from 
N(> F) toN(> Fexp[f + 5t]) ~ N(> F)[l-a(f + 6r)]. 
Here, a = —dlnN/d\nF is the (negative) slope of the in- 
trinsic galaxy luminosity function N(> F) and we have 
assumed r <C 1. Thus dust inhomogeneities induce a 
fractional fluctuation — cxSt in the galaxy number den- 
sity. Since St is correlated with the matter density field, 
dust extinction induces a correlation between foreground 
galaxies and background galaxies (quasars) such that 
w fb(0) = -a(ST{9')5 f g (6' + 6)). Here, S f g is the fore- 
ground galaxy number overdensity. On the other hand, 
lensing in duced fluctuations in galaxy num ber density are 
2(a — 1)k (Ba rtelmann fc Schnei der 2001), where the — 1 
term accounts for the fact that lensing magnifies the sur- 
face area and thus decreases the number density. Because 
of the different dependence on the slope a the signal of 
extinction and lensing can be separated simultaneously. 

The SPSS galaxy-quasar cross correlation measure- 
ment (Scranto n et al. | [2005) is consistent with the a — 1 
scaling and thus the dust contamination, if any, remains 
sub-dominant. Our dust models are consistent with this 
measurement, since the expected fractional contribution 
from dust extinction is — a/[2(a — 1)](StS^)/{kS0 ~ 
-a/ (a - 1) x (0.27,0.11,0.08,0.03) for dust model A, 
B, C and D respectively. However, such measurement 
is already at the edge of providing interesting dust 
constraints. For instance, model A induces at 8 = 
0.01° a negative correlation with amplitude ~ 0.003a6 g , 
where b g is the SDSS galaxy bias. This signal is al- 
ready very close to the measurement uncertainty (Fig. 7, 
IScranton et al. I (|2005h and the averaged (a) ~ 1 from 
their table 2). In principle, by combining color and 
flux dependences of the galaxy-quasar cross correlation 
and the color-galaxy cross correlation, it will be pos- 
sible to separate the contribution of lensing magnifica- 
tion, gray and reddening simultaneously 2 . The next 
generation of galaxy surveys such as LSST, ALPACA 
or PanSTARRS will provide foreground galaxy-quasars 
measurements that can achieve a S/N ^> 10. This will 
allow to discriminate the above dust models unambigu- 
ously, thus providing accurate constraints on the cosmic 
dust extinction and clustering properties. 

4. CONCLUSIONS 

Several new methods have been proposed for inferring 
the lensing magnification signal from a variety of cor- 
relation measurements. These involve SN la flux, the 

2 Private communication with Brice Menard. 
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fundamental plane and Tully-Fisher relation of optical 
galaxies. In this paper we have shown that contamina- 
tion of cosmic dust extinction may severely degrade such 
measurements. As an example inhomogeneities in the 
cosmic dust distribution may limit the S/N of SN lens- 
ing measurements to < 10 level. 

Billions of galaxies can be detected/resolved by the 
Square Kilometer Array 3 through the 21cm hypcrfinc 
transition line emission which is not affected by dust ex- 
tinction. In such a case the only scatters other than 
intrinsic ones in the Tully-Fisher relation (L ex vf) are 
induced by lensing magnification, L — » L(1+2k). There- 
fore lcnsing reconstruction using these galaxies is an at- 
tractive possibility, since it is free of some systematics 
associated with cosmic shear such as shape distortion in- 

3 SKA: http://www.skatelescope.org/ 



duced by point spread function. 

On the other hand, spatial correlation measurements 
of SN fluxes or galaxy-quasar will constrain the amount 
of cosmic gray dust and its clustering properties to high 
accuracy. This will not only provide a better understand- 
ing of IGM dust physics, but also a valuable handing of 
dust contamination in the SN la Hubble diagram. 
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